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Abstract. The northeast model is a spin system on the two-dimensional 
integer lattice that evolves according to the following rule: Whenever a 
site's southerly and westerly nearest neighbors have spin 1, it may reset 
its own spin by tossing a p— coin; at all other times, its spin remains 
frozen. It is proved that the northeast model has a phase transition at 
p c — 1 — p c , where /3 C is the critical parameter for oriented percolation. 
For p < p c , the trivial measure 80 that puts mass one on the configura- 
tion with all spins set at is the unique ergodic, translation invariant, 
stationary measure. For p > p c , the product Bernoulli-p measure on 
configuration space is the unique nontrivial, ergodic, translation invari- 
ant, stationary measure for the system, and it is mixing. For p > 2/3 it 
is shown that there is exponential decay of correlations. 

Key words: northeast model, facilitated spin-flip system, oriented perco- 
lation, exponential mixing 

1. Introduction 

The northeast model is the simplest nontrivial facilitated spin-flip system 
on the two-dimensional integer lattice Z 2 . Interest in such systems dates 
to jSj; for a recent review see [3]. The northeast model is specified by the 
following rules: Spins take values in the two-element set {0, 1}; the spin at 
site (x, y) € 1? may flip only at times when the spins at sites (x — l,y) and 
(x, y — 1) are both 1; and at such times the flip rates are 

(1) — > 1 at rate p; 

1 — > at rate 1 — p. 

A rigorous construction of the process is outlined in sec. |21 It will also be 
shown that the product Bernoulli-p measure u p on the configuration space 
is invariant and reversible. Note that unlike many growth models (e.g., 
the contact process) the northeast model is not additive (in the sense of 
Liggett Definition III. 6.1) or even monotone (in the sense that stochastic 
monotonicity is preserved by the evolution). 

The northeast model is the natural two-dimensional analogue of the one- 
dimensional east model studied by Aldous and Diaconis pQ. The east model 
is the spin-flip model on the integer lattice Z in which the spin at site x may 
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flip only at times when the spin at site x — 1 is 1 ; the flip rates at such times 
are given by ([T]). The main results of are rigorous bounds on relaxation 
rates to the product Bernoulli-p measure that are sharp in the limit p — > 0. 
The northeast model differs fundamentally from the east model in that the 
product Bernoulli-p measure v p is not even ergodic for p < 1 — (3 C , where (3 C 
is the critical parameter for oriented site percolation. The reason for this 
is apparent: If p < 1 — (5 C then v„— almost every spin configuration must 
contain infinite southwest clusters of 0s; these clusters must remain frozen 
at spin forever. 

We shall establish in this paper that the northeast model exhibits a phase 
transition at p c := 1 — j3 c by showing that for p > p c the product Bernoulli-p 
measure is ergodic and mixing: 

Theorem 1. The only ergodic, translation invariant, stationary measures 
for the northeast model with flip parameter p > l—f3 c are u p and uq. The only 
ergodic, translation invariant, stationary measure for the northeast model 
with flip parameter p < 1 — (3 C is vq. For p > p c , the product Bernoulli 
measure v p is not only ergodic but also mixing. 

Here mixing means that correlations decay to zero as t — > oo. Following is 
a more precise statement of what we shall prove. For any spin configuration 

denote by £ A its restriction to the set AcZ 2 . Similarly, denote by f A the 
product Bernoulli-p measure on the restricted configuration space {0, 1} A . 
Denote by £t the (random) spin configuration at time t in a realization of 
the northeast model. 

Theorem 2. If p > 1 — f3 c then v v is mixing for the northeast model in the 
following sense: For u p — almost every initial configuration a, and for every 
finite set A C Z 2 and every configuration £ A G {0, 

(2) iim p;fe A = c A } = ^ A (c A ). 

In light of the Aldous-Diaconis results, it is natural to ask if correlations 
decay exponentially for p > p c . We conjecture that they do. We have been 
able to prove this only for p > 2/3. 

Theorem 3. For each p > 2/3 and each finite set A of sites there exist 
constants C\,a^ > such that for any two configurations C,^ } n^, Vi > 0, 

(3) \P^{it = CVo = V A } ~ ^(C A )Mn A )\ < C A exp{-a A t}. 

We conjecture further that the exponential decay parameters a A = a A (p) 
can be chosen so as not to depend on the sets A. If this is true then for p > p c 
the northeast model has a positive spectral gap (see P), and the nature of 
the phase transition would be reflected in the behavior of the spectral gap 
as p ->■ p c . 

Sections [2HU are devoted to the proofs of Theorems Theorems ^ and 
12 will be proved in sees. [2HS1 and Theorem 01 in sec.E] In the final section 
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El we state a conjecture about the propagation of influence in the northeast 
model. 

2. Construction 

Two constructions of the northeast model are possible, one using the 
Hille-Yosida theorem (see Liggett jl], ch. 1), the second using a system of 
independent marked Poisson processes. The first has the advantage that 
it gives simple characterizations of invariant and reversible measures. The 
second yields detailed information about the time evolution of the process. 

2.1. Generator and Semigroup. Let X be the space of spin configura- 
tions on the lattice Z 2 , and for any subset A C Z 2 let X A be the set of 
configurations on A. Denote by C(X) the space of continuous real- valued 
functions on X (relative to the product topology), and by C*(X) the subset 
consisting of those functions that depend only on finitely many coordinates. 
Note that C*(X) may be naturally identified with U\C(X ), where the 
union is over all finite subsets A of Z2. For any configuration a, define A(a) 
to be the set of sites that are flip-eligible in configuration a, that is, those 
sites x whose nearest neighbors to the south and west both have spin 1 in 
configuration a. Also, for any site x, denote by <7+ and a~ the configura- 
tions that agree with a at all sites y 7^ x and have values 1,0, respectively, 
at site x. 

For any function / G C(X) and each site x € Z2 define 

(4) V x f(a) = (/(*+) - f(a~))p if a = a~; 

= (/(O "/(<£))(! -P) if a = <7+ 

For any / G C*(AT) define 

(5) Cf(a) = V a ./(a)l{x e A(a)}; 

Observe that because / € C*(X) depends on only finitely many coordinates 
the sum is finite, and hence the definition is a valid one. It follows from 
jlj, sec. 1.3 that the operator C on C*(X) extends uniquely to a Markov 
generator, and that C*(X) is a core for the generator. (Note: [1] uses a 
larger core, but the graph of C restricted to C*(X) is easily seen to have 
closure equal to that of C restricted to the larger core of 4 .) The Hille- 
Yosida theorem therefore implies that the closure of C generates a unique 
Feller Markov semigroup on C(X). This implies the existence of a Markov 
process with flip rates ([T|). 

2.2. Stationary and Reversible Measures. A Borel probability measure 
[i on X is stationary for the Markov semigroup generated by C if and only 
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if for every / G C*(X), 

(6) fcfdfi = 0. 

A stationary measure fj, is reversible if and only if for all pairs f,g£ C*(X), 

(7) j fCgdf, = J gCfdfi. 

These criteria yield easy proofs that the product Bernoulli-p measure v p is 
stationary and reversible for the northeast model with flip parameter p. By 
© and (|B|). to prove that v p is stationary it suffices to show that for any 
function / : X — > 1R that depends only on the spins in A and each x 6 A, 

/ V x f(a) du p (a) = 

J{aeF(x)} 

where F(x) = {a : x E A(a)} is the set of configurations a for which site 
x is flip-eligible. The event {x € A(a)} = {a € -^(x)} does not depend on 
the spin at x, and so under v p this event is independent of a(x). Moreover, 
the event a = occurs if and only if cr(x) = 1, and hence is independent 
of {a £ F(x)}. Since a(x) = 1 with probability p under u p , Q) implies that 
(with q = 1 — p) 

V x f(a) dv p (a) 

F(x) 

= [ ((/(<£) " ffa))pq + (f( a x) ~ f(<r+))gp) du p (a) 

Jf(x) 

= 0. 

Similarly, to prove that v v is reversible, it suffices to show that for any pair 
of functions /, g £ C{X) that only depend on the spins in a finite set A, 

/ g{oWxf{o)dv p (o) = I f{a)V x g{a)dv p {<j). 

JF{x) JF{x) 

Since the spin at x is independent of the event a G F(x), this identity follows 
by a simple calculation similar to that above. 



2.3. Construction via Marked Poisson Processes. Let P) be a 

probability space on which are defined countably many independent rate-one 
Poisson processes and countably many independent uniform-(0,l) random 
variables. Assign one Poisson process to each site i£Z 2 , and one uniform 
r.v. to each occurrence in each Poisson process. The occurrences in the 
Poisson process attached to site x mark the times of flip opportunities for 
site x: at each such time, site x queries its neighbors to the south and west 
about their current spins, and if both of these are 1 then x resets its spin 
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according to the value U of the uniform r.v. attached to the occurrence. 
The reset rule is 

(8) U < p => reset spin to 1; 
U > p ==> reset spin to 0. 

Note that such a reset does not necessarily change the spin at x: if the reset 
occurs at a time when the spin at x is (respectively, 1), then the chance 
that the spin changes is p (respectively, q). Thus, the flip rates at site x will 
agree with the specification (JX|). 

In order that this construction uniquely specify the spin at site x at each 
finite time, it must be the case that the spins of the southwest neighbors of x 
at query times can be determined. For this to be the case, the initial config- 
uration £0 must be specified. Throughout the remainder of this section and 
section [231 we assume that £0 = C f° r some fixed, nonrandom configuration 
(. (We suppress the dependence of the process £t on £, but the reader should 
note that the algorithm specified below will give a different time evolution 
£t for every different choice of £.) To determine the spin at site y at time 
t, one must search backward in time for occurrences in the Poisson process 
attached to site y: If there are no occurrences in this process between times 
and t, then the spin at site y at time t is set to its value at time 0. Other- 
wise, at the last occurrence time s before t, the backward search moves to 
the south and west neighbors y', y" of y and proceeds recursively. If it can 
then be determined that the spins at y' and y" at time s are both 1 then 
the spin at y at time t is determined by the attached uniform according 
to the rule (JSJ. Otherwise, the query continues to the last occurrence in 
the Poisson process at y before time s. To show that this algorithm ter- 
minates, we must show that the backward tree of queries initiated at (y, t) 
is almost surely finite. But this follows because the backward query tree is 
stochastically dominated by a simple binary fission process with fission rate 
1. (At each time a query is made, two new query processes are engendered. 
The additional queries in these offspring processes are mutualy independent, 
except when they coincide by virtue of a merger, e.g., SW and WS). 

2.4. Regeneration at Reset Times. Assume henceforth that the north- 
east process £t has been constructed as in sec. I2.3| using auxiliary marked 
Poisson processes and a specified initial configuration £o- For each site x, 
define 

(9) t x = time of first reset at x. 

Note that r x may take the value +00. Note also that r x need not be (and 
generally will not be) measurable with respect to the a— algebra generated 
by {&}o<t< 

Let < and < be the natural weak and strong partial orders on Z 2 : that 
is, x < y (respectively, x < y) if each coordinate of x is < (respectively, <) 
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the corresponding coordinate of y. For each site x, define 

(10) Ax = {y : y < x} and 
Bx = {y ■ x % y}. 

Sites in A x are those that may influence the times at which spin resets at 
x may occur; and sites in B x are those whose reset times are not influenced 
by site x. Let Nf be the the number of occurrences in the Poisson process 
attached to site x up to time t, and let {Uj}j>\ be the uniform-(0, 1) ran- 
dom variables attached to the occurrences of this Poisson process, listed in 
chronological order. Define a— algebras 

(11) Qt = *(N?) a < t ; 

Uf = <j(Uf)i< N x; 

If = a(U yeAx J^); 
T? = a(U x l)^ 
J t x = a(U yeBx Pf); and 

K? = *{j*ug*). 

Note that If° is generated by everything that happens in the entire con- 
struction up to time t. 

Proposition 1. For each site x and each time t > 0, the conditional dis- 
tribution of the spin £t(x) given the a— algebra JC^ on the event {t x < t} is 
Bernoulli-p, that is, 

(12) P{l- t (x) = l\K* Q )=p on{r x <t}. 

Proof. Observe that the reset time t x , and in fact the entire string of flip 
opportunity times at site x up to time t, are measurable with respect to the 
a— algebra Kf . On the event t x < t, there is a last flip opportunity at x 
before time t, and this occurs at the kth occurrence of for some k > 1. 
The spin at x at time t is determined by the uniform r.v. U%. But this 
uniform is independent of /C^. The assertion (|12|) follows. □ 

Corollary 1. If, for every site x in a finite set A, 

(13) P{t x < oo} = 1 

then the joint distribution of the restricted configuration & converges weakly 
to the product Bernoulli-p measure v£ as t — > oo. 

Proof. Consider first for illustration the case A = {x, y} where x is the 
westerly nearest neighbor of y. Observe that /Cf C 1C\ , and recall that 
it{x) is measurable with respect to /Cf . By hypothesis, P{t x < t} w 1 and 
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P{t v < i} ~ 1 for large t. Consequently, by Proposition ^ (first for y, then 
x), for large t, 

P(Zt(x) = Uv) = 1) = EE(P^ t (y) = 1 1 | Kf ) 

* P EP(Z t (x) = l\K%) 
«p 2 . 

The general case is proved by induction on the size of A. Fix a site y € A 
so that A \ {y} C B y : since A is finite, it must contain such a site. Set 
A' = A \ {y}. Observe that for each x E A', the spin £t(a?) is measurable 
with respect to 1C\. Let £ A be any spin configuration on A, and let C R be 
its restriction to a subset R of A. Then by Proposition ^ and the induction 
hypothesis, for large t, 

P(Zt = C A ) = EP&iy) = C y I )lfe A ' = C A '} 

« ^(CW^' = C A '} 

«^(C)^'(C A ') 

= ^(C A ). 

□ 

3. Stationary Measures, Ergodicity, and Mixing 

3.1. Examples of Stationary Measures. We have shown that the prod- 
uct Bernoulli-p measure v v is stationary for the northeast process. This is 
not, however, the only stationary measure: for instance, the measure vq that 
assigns probability one to the configuration with all spins is also station- 
ary, and ergo die. In fact, there are infinitely many distinct (and mutually 
singular) ergodic stationary distributions. A denumerable family may be 
built as follows. 

Let r be an infinite, connected subset of Z 2 such that 

(a) for every site y £ T, at least one of the southerly or westerly nearest 
neighbors of y is also an element of T; and 

(b) every connected component of Z 2 \ F is finite. 

Consider the initial configuration £o in which every site in T is assigned 
spin and every other site spin 1. If the northeast process started in this 
configuration then every site in T will remain frozen at spin 0, because any 
query tree (see sec. 12.3(1 that begins at such a site will have leaves in T. 
Similarly, every site x one of whose southerly or westerly nearest neighbors 
is in r will remain frozen at spin 1 forever. Now consider the evolution £^ 
in a connected component A of Z 2 \ F: because A is finite, and because its 
southwest border consists of sites that must remain frozen at spin 1, the 
process £ A is an ergodic finite-state Markov process, and so the distribution 
of £ A must converge as t — > oo to a stationary distribution /i A . The product 
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Ar of these measures // over all components A with the point mass on the 
zero configuration in T is a stationary distribution for the northeast process. 
Since each of the factors /U A is stationary and ergodic for the restricted 
process £^ , the product measure Ar will be ergodic for the northeast process. 
Observe that the measures Ar and Ar' will be mutually singular if T ^ V , 
because different sets of sites (T and V) remain frozen forever. 

The preceeding construction shows that there are infinitely many ergodic 
stationary distributions for the northeast model. None of these distributions 
is translation invariant. It is natural to ask if there exist translation invariant 
stationary distributions other than the product Bernoulli measures v v and 
v§. The answer is yes, as the following argument shows. Fix an integer 
m > 2, and let Tq be the set of sites with at least one coordinate is divisible 
by m. For each x = (xx,Xq) € 1? with coordinates satisfying < X{ < m, 
let T x = Tq + x be the translation of Tq by x. Note that each of the sets T x 
is invariant under translations by (m, 0) and (0, m), and that T x + y = T x+y 
where the addition in the subscript is done mod m. Define 

( 14 ) ^=^E A r- 

Then \i is certainly translation invariant, because a translate of T x is another 
r^, and n is stationary for the northeast process, because it is a mixture of 
stationary distributions. Note that m is not ergodic, because the initially 
frozen set T x must remain frozen forever. 

3.2. Characterization of Ergodic, Translation Invariant Measures. 

We now show that if \x is an ergodic, translation invariant stationary dis- 
tribution for the northeast model with flip parameter p then either \x = v p 
or fi = vq. Denote by P^ = P x fx the probability measure on 17 x {0, l} z2 
according to which the initial configuration £o is chosen randomly from \i 
and the marked Poisson processes used in the construction of sec. 12.31 are 
built on J 7 , P). Under P„, the distribution of £o is translation invariant, 
and so there exists a constant < r < 1 such that Pn{^o(x) = 1} = r and 
hence (by stationarity) 

(15) Ppittix) = 1} = r VxGZ 2 ,Vt>0. 

If r = then fi = vq. If r = 1 then at every rational time, every site is in 
spin 1; but this is impossible unless p = 1, because it would imply that no 
site ever flips to spin 0, as spin values are held for time intervals of positive 
duration. Thus, if < p < 1 then either r = 0or0<r< 1. 

Suppose that < r < 1. Since \i is assumed to be ergodic, the Birkhoff 
ergodic theorem implies that the long-time average spin value at any site 
x converges to r, almost surely. Thus, with P^— probability 1, site x flips 
its spin infinitely often, and so Pu{t x < ooVx} = 1. Hence, Corollary ^ 
implies that for \i— almost every initial configuration £ the conditional joint 
distribution of any finite spin block £^ given = C converges weakly as 
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t — > oo to the product Bernoulli measure Vp . It now follows that r = p and 
H = v p . 

3.3. Ergodicity and Mixing of v v for p > 1 — (3 C . To complete the proof 
of Theorem ^ and the first assertion of Theorem [2] we must show that if 
p > 1 — P c then Up is ergodic and mixing for the northeast model. The key 
is that if p > 1 — P c then for v p — almost every spin configuration £ there 
are no infinite southwest clusters of Os. In particular, the size of the cluster 
containing the origin is almost surely finite. 

Lemma 1. Define Mf to be the number of flip opportunities at site x up to 
time t, and denote by xs, x\y the southerly and westerly nearest neighbors 
of x. Then for every initial configuration 

(16) P(M X (oo) < oo and M Xs (oo) = M Xw (oo) = oo | £ = C) = 0. 

Proof. It is enough to show that almost surely on the event 

G := {M xs (oo) = M Xw (oo) = oo} 

the spins at sites x$ and xw will both be 1 at arbitrarily large times t. 
Denote by Ty the time of the kth flip opportunity at site y. On the event 
G it must be that Ty < oo for all k > 1 and both y = xs,xw- At each 
time Ty, the spin at site y is reset by a toss of a p— coin. Since the spin 
values at y = xs and y = x\y play no role in determining the reset times 
at either xs or xw, it follows that the pair of spins at xs,xw will be reset 
by independent p— coin tosses infinitely many times. Hence, the spins will 
both be 1 at indefinitely large times, w.p.l. □ 

Assume now that the initial configuration £o is chosen at random accord- 
ing to the product Bernoulli-p measure v p , independently of the marked 
Poisson processes used to determine the time evolution, and assume that 
p > 1 — C . To show that the process £t is mixing (and therefore also er- 
godic), it suffices, by Corollary Q to show that t x < oo almost surely for 
each site x. We will show that in fact M x (oo) = oo a.s. 

Suppose to the contrary that for some site x there is positive probability, 
say g > 0, that M x (oo) < oo. Then by Lemma El for at least one of the 
sites y = xs,xw it must be the case that M y (oo) < oo. It then follows 
that for at least one of these two sites, the spin at this site stabilizes at 0, 
because if £t(y) = 1 eventually and the other southwesterly nearest neighbor 
y' has infinitely many flip opportunities, then S,t(y') = 1 for arbitrarily large 
times t and so x must have infinitely many flip opportunities, contradicting 
the hypothesis that M x (oo) < oo. By induction, it follows that there is an 
infinite sequence y n of sites such that (i) each y n +i is either the southerly 
or westerly nearest neighbor of y n ; and (ii) the spin at each y n eventually 
stabilizes at 0. 

This contradicts the hypothesis that p > 1 — f3 c : In particular, it implies 
(by translation invariance) that for each K < oo there is probability at least 
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q that at large times t the origin will belong to a size-> K southwest cluster 
of sites all with spin 0. Since by stationarity of v v the size of the southwest 
cluster of 0s containing the origin has the same distribution at time t as at 
time 0, this is impossible. □ 

4. Exponential Decay of Correlations 

In this section we shall establish the exponential decay of correlations © 
for all parameter values p > 1/2. Without loss of generality, we may assume 
that the finite set A in © is a square, since every finite set is contained in 
a square; and furthermore, we may assume that the northwest corner of the 
square A is the origin (0,0). 

Denote by <9_A the exterior southwest boundary of A, that is, the set of 
all sites not in A that border A to the south or west. At any (Markov) time 
when the spins at all sites in <9_A are 1 it is possible for the sites in A to 
begin to reset, starting from the southwest corner and proceeding north and 
east. In particular, it is possible that all spins in A will flip to 1 before any 
of the spins in d-A reset. Moreover, at any (Markov) time when all of the 
spins in A are 1, it is possible that the spins will then reset one at a time 
starting from the northeast corner and proceeding to the southwest, in the 
order indicated below (for a 3 x 3 square): 

4 2 1 
7 5 3 
9 8 6 

Call such an occurrence a (total) A— reset. Define 

(17) ta = time of first A-reset. 

Observe that at time ta the spin at the southwest corner of A has just 
reset. At this time all spins in A have been reset at least once, in order 
starting from the northeast corner, and so the conditional distribution of 
the configuration £^ , given the evolution to time ta of the process in the 
region southwest of A, is the Bernoulli product measure v^. 

Lemma 2. To prove J2J) it suffices to prove that Vi > 0, 

(18) Pp p {r A >t}<C A exp{-a A i}. 

Proof. This follows from Proposition ^ by a similar argument as in the proof 
of Corollary^ □ 

Whether or not a A— reset will occur quickly depends on the current values 
of the spins at the sites to the southwest of A, in particular, on those that 
directly border on A. Since any such site with spin value will freeze the 
spins at its northern and eastern neighbors, it is impossible for A to reset 
until all the southwest boundary sites have assumed spin value 1, unless 
all spins in A are initially at spin 1. Similarly, any of the boundary sites 
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initially at spin cannot flip until their southwest neighbors have assumed 
spin value 1, and so on. Thus, at any time t the obstruction to a A— reset 
is the union of the oriented 0— clusters attached to the southwest boundary 
of A. Denote this union by K^: that is, is the set of all sites y for 
which there exists a (northwest oriented) path of 0s in the configuration £ t 
beginning at y and ending at the southwest boundary of A. Fix > 0, and 
define stopping times 

T < S < Ti < Si < • • • 

by setting 

T = 0; 

S n = mm{t > T n : K t A = 0}; 
T n = min{t > S n _i : K A ^ 0} A U, 
where t* is a fixed, nonrandom time. 

Lemma 3. Assume that p > 2/3. Then there exist constants C,a > 0, 
possibly depending on the size of the square A and on the choice of the 
constant t*, such that for all n and all t > 0, 

(19) P(S n -T n >t\l%)< Cexp{-at}. 

Proof. Let X t := |-fT t A | denote the cardinality of the 0— cluster set attached 
to A at time t; then S n is the first time after T n when X t = 0. At any time 
when the cluster K A is nonempty, its exterior southwest boundary consists 
of B > 1 sites, all at spin 1, and its interior southwest boundary (the set of 
sites that are flip-eligible) consists of A > 1 sites, all at spin 0; since each 
interior boundary site is bordered by at most two exterior boundary sites, 

(20) B < 2A. 

The cluster can grow or shrink by only one site at a time, either by addition 
of a site in the exterior southwest boundary or deletion of a site in the interior 
southwest boundary. Since there is always at least one interior boundary 
site that is flip-eligible, the rate at which jumps in X t occur is at least p. 
Moreover, at any addition/deletion event, the conditional probability that 
the event is a deletion, given the current configuration, is at least 

pA/(pA + (1 -p)B) : 

this is because the instantaneous rate of deletions is pA, whereas the in- 
stantaneous rate of additions is at most (1 — p)B. (All of the spins in the 
exterior southwest boundary of K A are 1, but some of these may be frozen.) 
Since B < 2A and p > 2/3, 

P A/(pA+ (1 -p)B) > 1/2. 

Thus, the jump process X t is stochastically dominated by a reflecting ran- 
dom walk on the nonnegative integers with negative drift. It follows rou- 
tinely that the first passage time to has an exponentially decaying tail. □ 
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Proof of j .?#)) ■ Let G n be the event that there is a A— reset between times S n 
and T n+ i. Since at time S n the 0— cluster set Kg is empty, the conditional 
probability of G n , given the history of the process £t up to time S n , is 
bounded below by a constant 5 > 0. (Recall that the spins in cLA must 
stay fixed at 1 until there is an occurrence in one of the Poisson processes Nf 
attached to a site y € cLA, regardless of what happens elsewhere.) Hence, 
by Lemma |3 

H 

(21) TA<S + Y, R i 

3=1 

where So, H, and Ri, R2, ■ ■ ■ are mutually independent, H is the number of 
tosses of a 5— coin until the first Head, and 

P{Rj >t} = Cexp{-at} V t > logC/a. 

Elementary arguments show that the random variable on the right side of 
(|21j) has an exponentially decaying tail. □ 

5. Influence Propagation 

A key component of the Aldous-Diaconis argument (although not explic- 
itly noted as such) is a proof that the influence of a spin 1 is propagated 
forward at a definite linear rate. In two dimensions, propagation of influ- 
ence is a much more subtle business, as there is no longer a unique path 
along which it occurs. A natural way to measure influence propagation in 
the northeast model is by use of influence regions Rt defined as follows. 
Consider an initial configuration in which all sites in the first quadrant are 
initially set to spin 0, and all other sites are initially set to random (by 
independent p— coin tosses). For each time t > define 

(22) Rt = {sites in quadrant 1 that are flipped at least once by time t}. 

Figure 1 below shows the results of a simulation (done by Marc Coram 
of the University of Chicago) in which sites in quadrants 2, 3, and 4 are 
initially determined by independent p— coin tosses, and p = .8. The figure 
shows the state of the system at time t = 1044.65; the size of the box shown 
is 500 x 500. The influence region Rt is shown in white, and the blue region 
Qt consists of those sites that have been "queried" at least once by time 
t. See section |21 for an explanation of the query process. The black specks 
indicate those sites in Rt or Qt whose spin values are at time t. 

The results of this and other simulations suggest the following conjecture. 

Conjecture 1. Forp > p c , the influence region Rt as a definite, nonrandom 
limit shape S. That is, with R denoting the union of all unit squares centered 
at points of R, 



(23) 
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Figure 1. Influence Propagation in Northeast Model 

in the Hausdorff metric, for some nonrandom set S contained in the first 
quadrant ofM 2 . 

Note that because the Northeast proces is not additive (in the sense of 
the usual methods for proving shape theorems (see, e.g., [2] and [H]) are 
not applicable. 
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